1907MNRAS..67..4271 


May 1907. Computation of Secular Perturbations. 


427 


Computation of Secular Perturbations. By K. T. A. Innes. 


I. 


It is well known that the Lagrangian equations for the variation 
of a planet’s elements (the planet being considered of infinitesimal 
dimensions) are exact. It is in the integration of these equations 
that difficulties arise. To avoid these difficulties several assumptions 
are made, such as that the planets can be considered in pairs (i.e. the 
problem of three bodies) and do not react on each other; that the 
elements on the right-hand side of the equations are constant, and 
that the perturbative function can be developed in a converging 
series. When we are concerned with the mutual perturbations of the 
eight major planets, the above assumptions are sufficiently good; the 
errors introduced by the first two assumptions are eliminated later 
by taking into account the powers of the masses higher than the 
first, the third assumption is nearly justified by the smallness of 
the eccentricities and mutual inclinations of the eight major 
planets. In cases of difficulty, the Lagrangian equations are 
abandoned and the perturbations found by other processes; thus 
Hill found that the Lagrangian method as used by Le Yerrier for 
the theory of the motion of Jupiter and Saturn was insufficient, so 
that in his theory he adopted the processes developed by Hansen. 

Gauss showed that the secular part of the perturbations (to the 
first power of the masses) could be found without using a develop¬ 
ment of the perturbative function. The rationale of Gauss’s method 
is as follows:—Let us imagine that the perturbative function is 
developed in a series as follows, 

P 0 = T 1 + T s sin M + P c cos M + P 2s sin 2M + etc. 

where the coefficients P contain, besides certain quantities 
depending on the elements of the disturbed planet, the elements 
and position of the disturbing body. Multiplying each side by 
■c£M and integrating around the circumference gives 



Then writing P 1 =^ 1 +y? a sin M x +p c cos M 1 + etc., where the co¬ 
efficients p depend only on the elements of both planets and do not 
•contain the time explicitly, a second integration gives 



It is on the quantity p l that the secular perturbations depend. In 
the general case it is impossible to find this quantity by algebraical 
expansions. Gauss proved that one of these integrations could be 
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reduced to elliptic integrals; the second has to be a mechanical 
integration. Gauss, as is well-known, gave a geometrical title to 
his paper and did not prepare his formulae for numerical application. 
Bour and others have given further geometrical interpretations of 
Gauss’s formulae, which, interesting as they are, call for no attention 
here. 

Since Dr. G. W. Hill published his first paper on Gauss’s 
Method of Computing the Secular Perturbations of the Planets 
(Astron. Papers prepared for the use of the American Ephemeris, 
1882; Collected Works, T. ii., 1906) many numerical applications 
of this method have been made. Hill did not leave the method in 
final shape, for no sooner had he finished his paper than he pro¬ 
ceeded to give a substantial modification in an addendum by 
making use of Gauss’s arithmetico-geometrical mean. Both of 
Hill’s solutions require extensive tables of the elliptic integrals for 
their easy application. Hill provided these tables for his first 
solution, and Callandreau and the writer did the same for Hill’s 
second method. Previously to Hill’s work there had been but 
little use made of Gauss’s method. We learn from Bode’s Astro- 
nomisches Jahrbuch , 1819, p. 229, that before Gauss had actually 
published his memoir, Nicolai had computed the sec. var. of the 
elements of the Earth’s orbit by this method, but no details of the 
work were given. The late Professor J. C. Adams applied the 
method to the perturbation of the orbit of the November Meteors 
(see M. W., 1867, Apl.; Collected Works, T. ii. pp. 194-200). It is 
well known that the method as expounded by Gauss requires the 
solution of a cubic equation, but Adams avoids the cubic by neglect¬ 
ing a portion of the terms factored by the square of the eccentricity 
of the disturbing body, and so reduces the equation to a quadratic. 
This is no doubt the modification which Adams stated greatly 
facilitates the application of Gauss’s formulae. As a practical 
step this modification was quite justified, but it takes away from the 
interest of the problem, and, as will appear immediately, is of no 
real advantage. Halphen, in 1886 (C.R., T. ciii.) announced that 
the problem could be completely solved without the labour of obtain¬ 
ing the roots of the cubic ; he gave an exposition of his method in his 
Fondions Elliptiques , 1888, T. ii. I was unable to follow Halphen’s 
exposition, to which Tisserand’s explanation added nothing 
(Mec. Cel., T. i.), and it was only on obtaining by chance a copy 
of Dr Louis Arndt’s paper on Halphen’s solution (Recherches 
. . . . Pert. Seculaires, Neuchatel, 1896) that the solution was 
grasped. The formulae for numerical application given by Arndt 
assume that g z is always positive, so that his solution will only 
apply to the easier cases, although Halphen drew special attention 
to the case of g B being negative. Arndt has provided tables for 
the cases in which Legendre’s modular angle 0 will not exceed 45°. 
Hill has again revised his solution, and shows that by using Jacobi’s 
nome q tables are practically unnecessary, but the roots of the 
cubic are still required (Astf. Jour., No. 511, 1901). 

As already stated, Gauss’s method is rigorous to the first power 
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of the masses. It has been applied to most of the mutual actions 
of the major planets, but unnecessarily so in many cases, as the 
tables given by Newcomb and Le Yerrier are ample and much 
simpler to use unless the ratio of the mean distances exceeds one- 
half ; in other words, the small inclinations and eccentricities of the 
orbits of the major planets permit algebraical expansion, which 
only becomes long when this ratio exceeds one-half. 

This method has its r 61 e in computing the first order secular 
perturbations of the periodical comets and minor planets. An 
interesting application to Eros has recently been made by Herr 
W. Dziewulski ( SdJculare Marsstorungen . . . des Eros , Cracovie, 
1906). He shows incidentally that this method gives better 
results in the case of Eros than the Lagrangian method, which takes 
into account all powers of the masses, but not of the inclinations 
and eccentricities. In fact, the Lagrangian approximation will in 
some cases not even give the correct signs, far less the correct 
figures. Tolerable approximations are only given in the cases of 
the planets with preponderating masses. 

Mr C. J. Merfield, of Sydney, has recently announced the 
completion of a computation of the secular variation of Eros. 

In this paper Hill’s first exposition has been adopted as a 
working basis, and in the introduction of the elliptic functions 
Schwarz’s notation, as developed in Padd’s translation ( Formules et 
Propositions , Paris, 1894), has been followed, except that the 
quantity h is replaced by Jacobi’s letter q. No proofs of formulae 
found in previous papers are given. 

The aim throughout has been simplicity in numerical application. 
The final formulae have been framed for that end, and have been 
submitted to numerical tests. The pure mathematician dismisses 
one with a series guaranteed to be absolutely convergent for 
certain values of the argument, but a trial often proves that the 
series, though it may be absolutely convergent, is arithmetically 
absolutely useless. It will appear that if an accuracy of six places 
of figures suffices, the use of series can be entirely evaded, a result 
by no means to be expected a priori. 

The two main points in this paper are (1) the simplification of 
the formulae by the introduction of the Weierstrassian cubic; 
(2) the new series for the computation ,of the periods of the 
elliptic functions,—series which are remarkably inert, nearly all 
the variation being thrown on to a simple trigonometrical factor 
of the series. 

Thus the computer has the choice of two main methods in 
applying Gauss’s principle,—(1) solving the cubic and using 
Jacobi’s nome q (Hill, Ast. Jour., No. 511); or (2) using the 
invariants of the cubic and two hypergeometrical series. Either of 
these methods is much shorter than the earlier solutions; which of 
the two is the better in actual practice cannot be answered 
offhand. It may be remarked that the former uses almost entirely 
trigonometrical logarithms, which the latter does not, and this 
alone is a considerable advantage. 
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2 . 


All the equations and formulae wanted in practice are collected 
later on. Here it is unnecessary to write down the six Lagrangian 
equations; as a type we confine our attention to that for the 
eccentricity 

— = cos = - n CQS ^ Tsin yR + (cos v + cos E)sl 

dt Y dt i+ra L J 


We have to find that part of the values o£ R, S, and W which is 
independent of the time (or of the positions of the planets in 
their orbits). 

As already stated, two integrations are required, one of which 
will be done mechanically. The above equation is therefore 
written 



J ML mI cos <b sin E-R 0 - (~e - 2 cos E + - cos 2E N )?S 0 1 
i+m L r \2 2 Jr J 


The mechanical integration need only be made on the quantities 


-R„, -S 0 , and where 

r r r 


-R 0 = Af + — A^ c) cosE + i-Af sin E + — Af cos 2E + etc. 


with similar equations for —S 0 and -W 0 . 

r r 


Thus the final equations take the form 



1 _ m i n r 

1 

1 *+«* 
N 

_ 1 

J 00 i+mL 




The addendum to Hill’s first paper may be referred to for 
further matter concerning these equations. 

We have 



with corresponding expressions of the same form for 



and -W 0 . 
r 
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Putting 


a 
a n 


- = a 


7 

r 0 — a( 1 - e cos E) = — 


a, 


k cos (v + K) = A c 
- k sin (v + K) = B c 


sin n x 


= a 


k Y cos <j> l sin (/; + K x ) — A 
h\ cos (f> l cos ( v + K 1 ) = B, 
cos <j > 1 cos IIjl = C, 


A 0 =I+r o + 2e i r 0 A c 

B 0 cos € = r 0 A c + e x 
B 0 sin 6 = r 0 A s 


Note, Hill's A = a*A 0 , etc. 


= A 0 - 2B 0 cos (E a - e) + O 0 cos 2 E 2 

Hill's A = a Y A 0 


then 


-B 0 = — 

r 0 27r 


a‘ 


> A c (cos E x - + A s sinE x - r Q 


(1 - e x cos E 1 )^E 1 


-S 0 =~ 
r 0 27r 


0 B c (cos Ej - e x ) + Bg sin E t 


(1 - cos E^dEj 


/»27T 

1 . _ C,(cos E x - + C s sin E x Wlr 

-W 0 = — asm Jr 0 —- 1 -—^ s --(i-^cosE^E, 


If we remove the constant factors a 2 , a 2 , and a sin J respec¬ 
tively, the numerators of the fractions on the right may be written 
as 

/ _ g e L + [^( 1 + e[) -fe{\ cos E : + h sin Ej - he 1 sin E x cos E x - ge 1 cos 2 E x , 


where/, g , and h have the respective values 


o 

e 1 

! d 

4 -s 0 

a 2 r 0 0 

1 " W x 1 

a sin J r 0 0 r 0 

/ 

~ r o 

o 

o 

g 

A c 

Be 


h 

A s 

b 8 

a 
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3 - 

Altering the variable from E x to T changes the equations for 

-4 — R 0 , etc., to the form 
a 2 r 0 


i I\ + r 2 sin 2 T + r 3 cos 2 T 


dT 


27 r 


wherein 


is now G x - G 2 sin 2 T + G 3 cos 2 T 
' but we may now omit the subscript 0 as unnecessary. 


• • ' 7 T 

The elliptic integrals, using the limits ° to —, are 


dT 
A 3 ’ 


sin 2 TdT 

A 3 


and 


cos 2 TdT 
A 3 


(«) 


The G’s are the roots of the equation 

x 3 -(A- C)x 2 + (B 2 - AC)x + CB 2 sin 2 e = o 
or x 3 — ky + k 2 x — k 8 — o 
Putting 

= 0.-^1 


3 


e,= G,-A e.— - G,- 
2 2 3 3 


transforms 

Gj — G 2 sin 2 T + G s cos 2 T into e 1 - e 2 sin 2 T - e 3 cos 2 T 
Introducing the Weierstrassian functions by means of 

o, s = e 0 



s j n 2 rp _ 

. e i e z 

s - e 2 



e 2 e 3 

s - e 1 

gives 


cos 2 T — 

e 2 e i 

s-e s 

dnd 


e 2 ~ e 3 

s - e Y 


when T = < 7r 


- s = e a 


dT _ i /(e 1 - e s ) (e . 2 - e x ) i 

2 V (*-e 2 )(s-e 8 ) ‘ 

A 2 = («j — e 2 ) (ej — e 8 )/(s - 6j) 
dT i s - e. 


ds 


A* 2 (e, - e 2 ) (8, - e 3 ) ' V /( S - ^ (. - e 2 ) (. - e 8 ) 
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433 


Putting 


and noting that 


= 2 ^(pu - e x ) (J)w - e 2 ) (p« - e 8 ) 


we have 


— — ' 1 - du 

A 8 ( e1 — e 2 ) (e^ — 63) 

Having pu = e 5 when u — a/ and to e 2 when M = or w + a/ and 
writing JG = {e 1 - e 2 ) - e 3 ) (e 2 - e 3 ) and noting that 

, ,@1 

- pudu = a-^rU 


or 


so that 


we obtain 


j pudu— - 


pudu — r} 


©1 

© 


u 


7 r 

2 dT e 9 - e„ x 

— 3 ( ei o. + ,) 




sin 2 T(iT - e v , 

—rs— = -4=^(e 2 o> + 1 ?) 

A s 


(6) 


COS 2 T^T 6 9 - ^ . . 

A 3 + 

The proofs of the above formula will be found in Pad^’s or 
Schwarz’s paper already referred to. The italic G is not to be 
confounded with the three G roots. 

The last set of formulae correspond £o Dr. Hill’s cos 2 k 

m 6 

M 

and —g sin 2 k in the Amer. Journ. of Math., t. xxiii. pp. 321-2. 
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Should the cubic equation be solved, the numerical value of 
these formulae may be found as follows :— 



1 1 + 3 2 q 2 + 5 2 q 6 + 7 2 # 12 + etc - 

4 Jq (1 + q 2 + q Q + q 12 + etc.) 3 


Sq 


1 + 2 2 q B + 3 2 <? 8 4 - 4 2 g 15 + etc. 
(1 + 2 q + 2^ 4 + 2^ 9 + etc.) 3 



cos 2 TV 7 T _ Je 1 - e 2 
“A 3 “V 


Sq 


i - 2 2 q B 4 - 3 2 g 8 - 4 2 ^ 15 + etc. 
(1 - 2q + 2q 4 - 2q 9 + etc.) 3 


where q is to be derived as follows:—Let 


then 


——— = cos 2 0 and JZ 
e i e z 



2 


q = ¥ + 2 (¥) 5 + 1 5 W +15 o( JO 13 + efcc * 


or in practice 

log 2 = log + [9-9388] (iiy 
(see M.N., R.A.S., vol. lxii. p. 494). 

All these formulae are extraordinarily convergent. They can 
also be derived from the formulae given on p. 333 of Hill’s paper 
( A.J ., No. 511) in virtue of the relations 


T 1 + q 2 + q 6 + q 12 + etc. ~j 2 
— 4 S/ ( 1 [_ J + + 2^4 4. 2^9 4. et c> J 


i - 2q + 2<? 4 - 2 q 9 + etc.! 2 
_ 1 + 2q + 2</ 4 + 2 q 9 + etc. J 


2 
7r 


dT 


J1 - h 2 sin 2 T 


(1 + 2 q + 2 q A + 2 q 9 + etc.) 2 


7 r 

£ f 2 ^/i - & 2 sin 2 TO = 
”• J „ 


j _j_ ^ 2^2 -4- _f_ y2gl2 4. e tc. 

i - 3^ 2 + 5^ 6 - 7# 12 + etc. 


X 7^1 


the last formula being one of some half dozen representations of 
Legendre’s E by means of Jacobi’s q function. 
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Introducing the equations ( b ) into the equation (a) and putting 

M = (e 2 - e 3 )I\ + ( e 1 - e 3 )V 2 + (e 2 - e 1 )V 3 
N = efaz - e 3 )r i + e 2 (e ± - e 3 )T 2 + e 3 (e 2 - efT z 

enables us to write 

^ n ° = l^G [ Mj ? + Na> ] • • • ( c > 

with similar simple expressions for the two other components. If 
M and JST are written in determinant form they are 


M = 


Ti 

~ r 2 

-r. 

. N = 

e i r i 

~ e 2^2 

~ e 3^3 

e i 

e 2 

e 3 

e i 

e 2 

e 3 

I 

I 

I 


1 

I 

I 


5 - 

We now proceed to the elimination of the roots of the cubic 
(e 2 , e 2 and e 3 ) by means of the invariants g 2 and g 3 . The following 
formulae are collected here for easy reference:— 

+ G 2 - Gr 3 = &i. “ GjG 3 - G 2 G 3 = ft 2 . — 0 j 0 2 0 3 = Jc 3 

lg 2 = \ = k\- zK = - 3 («i e 2 + Vs + Vs) 
ids “ - 9 k A + 2 A) = W« 

i6G=gl - 2 jgl= i6(e x - e 2 ) - e 3 ) (e 2 - e 3 ) 

or for the last 




i, 1 




i> 1, 

I I 

e \ e 2 
e 2 e s e i e s 


1 

1 


e i e 2 


&i e 2 e 3 

We have now four determinants; let us multiply the first by the 
third, and then the second by the fourth. The first result is 
written at length, 


M = 


v -V -V 
A 1 A 2 8 1 

o 

3 


e 2 e 3 f l e 1 e 3 V 2 — e 1 e 2 T 8 , 
3 e i e 2 e s 

^e 2 + e 1 e 3 + e 2 e 3 


“ e 2-^2 “ e 3-^3 

— 2 + 6263 + e 2 e 3 ) 

o 


To abbreviate put 




r 2- 


v * ss 'l' 

e 2 e 2^l ~ e i e B^2 ~ e i e 2p3 = X • * ($) 

— * 2 r 2 — e 3 r 3 = </> 

and replace therein the constituents in e 1} e 2 , and e 3 by their 
appropriate equivalents in g 2 and g 3 ; then expanding we obtain 
finally 

M \JG = <f,yl + f(7 2 x - 1 -g t 4 > 

A similar process gives 

^ JC= — il/j^g 2 g 3 — f g s x + 
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w 

so that (c) becomes 

^ h(g 2 v- +3x) + - 99 s v)<)>] 

with similar equations for S 0 and W 0 . 


6 . 

Bruns first pointed out that the periods of the elliptic functions 
could be found without a knowledge of the three roots (H. Bruns, 
Perioden der Ellip. Integrate , Dorpat, 1875). 

This remarkable investigation does not deserve the neglect it 
has received at the hands of English mathematicians. Bruns does 
not adapt his formulae for numerical use, and I have found it 
better to adopt another absolute invariant. If we write 


r 08 

cos 1=3 Jl-j 

si 


(o°< t< 180 0 ) 


we have 


2 J_I/ 3 \ 51 1 -p/ 1 5 . ,t\ 

»■ JG 2 v- v 2 ^")“ 8 xi » F (e’ 6’ 2> sm 2 ) 

v COS 2 — 


2 
7r 






cos^ 


/ i 7 

— F - -7 2. sin 

i V 6 ’ 6 * 2 


2! 


Hence 


JL % _ I 

a 2 r 0“ 8 


cos^ — L 


^ + 3 x)s^’(|) + ^ 7 l'(- g)] 

with similar equations for S 0 and W 0 . The abbreviations fQ^ 
and F^- require no explanation. 


7 - 

It remains to find the values of i p, Xy aQ d <£• Replacing the 
e’s by the G’s in (d) we have. 


^»r x -r 2 -r 8 . 

* = ( G 1 ~ D ri ~ ( G * - 7 ) 1V * + ( G 3 + y) r 3 = G ! r i - G 2 r 2 + G 8 r 3 ~ 

X=- G 2 G,r 1 + G 1 G 3 r 2 - G 1 G 2 r 3 + - (AY ^ 

3 \3 / 
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£ 

e 1 s- 


d 
• ^ 

03 


z/i 

e I ^ 

w C 5 


pf 

«S I ^ 

i<M 

M I CJ 


U ^ 
I II 

rH 

P 


vu 

t d 

*03 

o 

p 

«o 

O 

+ 

a 

o 

+ 

o 


w 

d 

• rH 

02 

O 

PQ 

00 

p 

+ 

P 

p 

+ 

o 


w 
o O 

” *3 


o 


vu 

d 


* + 
O I 

1 + 


p 

co 

o 
+ ' 
rp *%■ 

o 


+ 

-e- 


P 

rH 

0 


H 

1 « 

o J~* 

•<3 "O' 

^ W 
O o 

V 0 

w u; 

.5 d 

50 'm 

P° W ° 

■ft 

of 


v2 

OJ O 

C 0 

CJ rH 

p -e- 

1 m 
1 o 
■ o 

u» 


<1 

CO 

P 

vw 


.a w 

CO O 

o %H 

p ® 

I 


<1 

CJ 

I 

O 

<1 

<5 


.a 

03 


P 


3 O 


02 

I 


P 

CM 

o 

rH 

O 


P 

CO^T-i _ 

I ON 


o 

+ 

pT 

10 

o 

Cl 

0 

I 


+ 

-e- 


co 
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£ 

o 


m 

d 

o 

• I™« 

•** 

o 3 

O 

Pi 

P 

fr. 

03 


c3 

o 

"E 

03 

s 

d 

d 

o 

£ 

■+* 

n 3 

d 

CJ 

J 

"d 

8 

t-i 

o 

*n 

&0 

d 

• rH 

o 


0 


+ 

03 

02 

0 

d 


d 

*—1 
CD 

• rH 

1 

C^> 

VU 

• rH 

a 2 

O 

O 

'P 

a 

0 

P^ 

03 

1-i 

£8 


r—H 

+ 

1-1 

d 

5 

Pi 


O 

II 

* 4-1 

p 

03 

03 

03 

^1 

03 

P 

03 

EH 
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9 - 

Let 

a = semi-axis major corrected for the constant part of its 
perturbations 
n = mean motion 
e = sin = eccentricity 
7T —long, of perihelion from equinox 
i — inclination of orbit to ecliptic 
12 = long. of ascending node 

X = long. of perihelion from a fixed point in orbit 

(X) = 7T — 12 

r = radius-vector 
E = excentric anomaly 
v — true anomaly 
m= mass (Sun’s mass unity) 

of the disturbed body; the addition of a suffix ( 1 ) indicates 
similar quantities appertaining to the disturbing body. 


10 . 

- sin sin J(<E> + \k) = sin 1(1^ - 12) sin + i) 

- sin \ J cos |(<t> + \k) = cos ^(Oj - 12) sin - i) 

cos | J sin |(d> - \k) = sin |(12 1 - 12) cos + i ) 

cos | J cos |(<t> - \k) — cos J(12 1 - 12) cos - i) 

check formulae 

cos p sin q = sin i 1 cos (12 t - 12) 
cos p cos q = cos i x 
cosp sin r = - cos ^ sin (12 x - 12) 
cosp cos v = cos (12 1 - 12) 

ship = - sin i x sin (Oj - 12) 

then 

sin J sin < 3 > = sinp 

sin J cos <t> = cos p sin (i — q) 

sin J sin — r) — sin p cos (i - q) 
sin J cos (q> — r) = sin (i — q) 
cos J = cosp cos (i — q) 

rr = oj - . u 1 = o) 1 -'& 

k cos (K - II) = cos II 1 . k x cos (Kj - II) = cos II 1 cos J 

k sin (K — II) = - sin cos J. k x sin (K x - II) = - sin IIj 

check ; k 2 + k^ = i + cos 2 J 

When the Earth is the disturbing planet, these equations simplify to 

J = i * 

<& = o 
II = 7 T — 12 

ITj — 7 Ti — 12 
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and when the Earth is the disturbed planet, to 

J = ^2 


#=i8o° 

IT = I 8o° + 7T — O x 
n x = i8o° + 7 r 1 - 0 1 


II. 


The various functions of E which follow should now be computed 
for 4j equal parts of the circumference. The larger the number of 
parts is, the more exact will be the results of the mechanical 
integration, but j should not be taken larger than necessary. A 
measure of the accuracy attained is given by adding the odd values 
and the even values of each of the quantities : if these differ, more 
values of E are required. Reference to Hill’s original example or 
to those of later computers will assist the judgment. 

Take a — a/a^, then omitting the suffix 0 used in the demon¬ 
stration, when no confusion will be created thereby, we have 

r 0 cos v = a(cos E - e ). r 0 sin v = a cos cj> sin E 


?’ 0 = a(i - e cos E). (check) 
k cos (v + K) = A c k x cos t/q sin {v + Kj) = A s 

- k sin (v T K) = B c \ cos ^ cos (r + 1 ^) = B s 


A=i+rJ+2e 1 r 0 A c . 


l'0 a « * 


B cos e = r 0 A c + e x . Bsine = r 0 A 5 

k 1 = A-e^. k 2 = B 2 - Ae*. - k 3 — ejB 2 sin 2 e 

|^ 2 = X = k\ - 3^ 2 g 3 = -gV(2^ - 9 kf 2 + 2 >jk 3 ) 


cost — , (o°<t<i8o°) 

cf 1 


cos 1 — 


D = (1 + e‘j)B cos e - <q(A + e J) 




+ B C D + B S B sin e 
+ sin n x D 4 - cos <^>2 cos n T B sin e 


= o 


= - sin e cos ^ cos J 


= - sin e cos <jfq cos (v + n) 
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With the aid of the tables or formulae compute 

|—— F(i/ 6 ,s/ 6 > 2 > Bin* i)-F A 

O o L \ 2 / 


and 


Then 


cos 2 — 


8 TT^ - 1/6 ’ 7/6 ’ 2 ’ sin2 2) = ^ B 


COS 2 - 


— r R 0 = f3X - 


a 2 r 
1 a 


A* 


a 2 r 

ii a 


-~S 0 = ( 3X ) » „ „ 


a sin J r 


- W 0 = [( 3X ) » n „]'o 


LXVII. 7, 


12 . 

With the 4/ values of these three functions, which may be 
designated 

R (0 >, R (1 >, R (2) , .... RW- 1 * 


S ,0 >, S (1 \. S w_1) 

. 

compute 


A^-i-jRW + RW + R® .... RW -«1 
4 / 

£ A‘ c) = — TR™ + R (1) cos £- - 4 - R <2) cos * - .... +RW- 1 >co 3 V ~ 1 E~\ 

2 1 47 L 72 .? 2 ^ 2 J 

£ A 1 * 1 = — f" R' 1 ' sin £ + RH, sin % . . . . + RW-D sin 4 ./~ £ 5 ] 

2 4 y l ^2 .? 2 ^ 2 j 

£ Aj’ = — Trio)+ RW cos i - + R <2) cos 4 - . . . . 4 -R w ~ 1) cos f Mz O il 

2 4iL 7 2 ^2 j 2} 

1 a£" = — T E' 11 sin % - 4 - R (2) sin 4 - .... 4 - R w-1) sin £Mz 2 )jf] 

2 4 ?'L J 2 ^2 / 2J 

with similar equations for B^ 1 .... and C^ 1 C[ c and C^. 

(Note :— Aj’ is not required.) 

For any special case of 7 these equations become greatly 
simplified : thus for j — 2 we have 

A ( ; ) = -J[R< 0, + R (1 > + R ,2) .... 4 -R' 7 '] 

JAf = -J[R (0) - R (4) 4 - (R (1) - R ,3) - R (5) 4 - R ,7 >) cos 45°] 
l Af = |[R< 2 ' - R<« 4 - (R (1 > 4 - R |3) - R' 6) - R' 7 ') sin 45°] 

\ Af = i[R 101 - R ,2) 4- R ,4 > - R |6> ] 

JA“ = i[V> - R (3) 4 - R ,4) - R <6 ’] 
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and so on. For other values of j the equations are easily formed, 
or may be found in Hansen’s Auseinandersetzung or elsewhere. 

The equation 

sin A[ s) + cos = o 
affords a useful check. 


13 - 


The final equations are 

m = «* T - (i + f) tf'-te 1 - Bf> + ~ i Bf] 

L dt Joo I+m l \ 2 / 0 2 1 2 2 2 J 

te] = 37 L a 2 r«A« cob * - I Af cos <t> + { 2 -e 2 )- Bf - i - 1 B<’>] 
LrfiJoo i+m 2 L 0 ^ 2 1 '2 1 2 2 2 J 

|te] = -- 1 -- a sin J cos <D |Y 2 Of - eC‘ c,V ) sec <f> - - Cf tan a,] 

LrftJoo 1 +m L\2 1 V 2 1 J 

tT—1 = « sin J cos co f—C! c| - eCf) sec <f> tan co + -Ic'* 1 "] 

L eft Joo 1 +m LV 2 V 2 1 A 

+ ^ +4 'T i " - ' 7 A “] 

Calculation of F^-i- y 2 , sin 2 -^ and F^ - -i-, 2 , sin 2 -^ 


sin % 


14. 


When l is less than 90° the calculation of these series is not 
long, hut if an accuracy of two units in the 7th place of the 
decimals of the logarithms is sufficient, the use of series can be 
avoided. This accuracy is indeed greater than either our know¬ 
ledge of the planetary masses or the precision of observations 
requires. Besides these,“the change in the secular perturbations 
due to higher powers of the masses than the first, is of a more 
important character. 

The approximate formulae are 
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442 


The two series are inert:— 


l 

F (t*1’ 2 ’ sin t) 

F ("T’i’ 2 ’ 

O 

O 

1*0000000 

I *0000000 

o 

90 

0*9415120 

1*04^2298 

180 0 

0'8i8siii=—-2 

r 2 0 

I’i 459 r 5 6 = — - z - 


IT 7 

* 5 


If l exceeds 90° the best procedure seems' to be as follows :— 
Take ^ = 180° —6, 

then proceeding at once to the actual quantities wanted and 
writing 

2 108 1 


Sill 

we have 


in t [ 7F (~ i’ i> 2 ’ sin t) + 5 008 I■ 2 ’ sin2 ir)_ 

5 2 , sin 2 i) =H- l log e l _L_ f(I , 2 , sin 2 h) 

cos 2 - V6 6 2> ^ cos 2 — V6 6 2> 

2 2 

-—F ( ~ 1 , 1 .2 , sin 2 = H cos 4 + - log, - - 1 — F ( -~, ? , 2 , sin 2 - 1 
,h 2 i V 6 6 2/ *■ 2„ ns 2h V 6 6 2 

To compute q put 


cos 2 - 
2 


COS 25 -i 
2 



,— tan -1 = cos 2y 
v 3 3 


and 

/s/tan y — cos 2 A 


then 

l = tan 2 A 


and 

q = 2 (^Z) 5 + etc. 

0 


In practice it may be assumed that log q = log \l. This 
assumption is in defect by the following quantities :— 

q Correction to log \l in units of the 7 th place of decimals, 
o° - 40° +o 

40° i 

5 °° 2 

6o° 4 

70° 9 

8o° 17 

90° ' + 3° 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ at NERL on June 28, 2015 




















1907MNRAS..67..4271 


May 1907. Computation of Secular Perturbations. 


443 


As lo£„ - 

oe q 

logarithms— 


~ is the natural log of ~, we have in ordinary 


log (L log e l'] = log of log - + 9-8650658 



15 - 

Numerical Examples. 


In the. Astron. Nadir., Iso. 4068, Dr Arthur B. Turner has 
computed the secular perturbations of Mars by Jupiter using 
Arndt’s formulae. 

The preceding formulae are applied to this case for the case 
of E = o°. We have J=i° 26', n=i49°*8, II 1 =i 88 0, 4, K and 
k i = 3 2I °'4, * and A i = I > log »•„ = 9-4241, v = o°, log A c = 9‘893, 
log B c = 9-795, log A s = 9’794«, log B 5 = 9-892, A 0 = 1 -09, 
log B 0 sin € = 9*219, kj — 1*088, k 2 = 0*090, * 3 = — o*ooo, g 2 — 1*218, 
^3 = °*2 5 i , l= 14,20, E(-J-) = 0*998, E( --J-) = i*ooi, D = 0*204, and 



4 > 

X 


0*002 


0*072 

0*061 


o‘i34 

0*050 


0*001 


With these figures Turner’s results are reproduced. 

The next example is taken from Herr W. Dziewulski’s paper 
on the secular perturbations of Eros by Mars. In the present 
position of the orbits for E=i20°, the podular angle 0=58°*9. 
In the notation of the present paper we have J=ii°*5, H = 6°*7, 
n i = 2i9°*6, K=i47°*6, K 1 =i46°*5, log* =9*996, log ^ = 9*995, 
log r 0 = 0*0268, v — i30°*6, log A c = 9*152, log B c = 9*992, 
log A s = 9990%, log B s = 9*084, A 0 = 2*i6, log B 0 sin €=o*oi6w, 
^1 = 2*151, *2=1*120, *3=-0*009, ^2 =i ‘ 6 9o ) y 3 =-0*300, 

t = 135°*!, t 2 = 44°* 9 , y=4o°*6, X = ii°*i, log JZ = 8*288, 



X °’ I 95 





leading to Dziewulski’s results. 


[Eor Tables, see paper by Mr Bobbins.] 


Johannesburg. 
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Mr Frank Bobbins , Tables for the 


LXvii. 7, 


Tables for the Application of Mr lnnes > s Method. 
By Frank Bobbins. 


As an appendix to the preceding paper on “The Computation 
of Secular Perturbations,” the author asked me to compute for 
each degree of the quadrant the logarithmic values (base io) of 
the two functions of iota required for the convenient application of 
his method. 

In the hypergeometric series F (a y x) in the first case 


a has the value - 

6 


e-i 


y — 2 


• o t 

# = sin J — 


and in the second case 


a = 


i 

6 



y — 2 



For convenience of designation the tables are headed Minus F 
and Plus F according to the sign of a. 

Vega’s (1794) ten-figure logarithms, corrected by collation with 
the copy in use at H.M. Nautical Almanac Office, were used, and the 
natural values of the individual terms were taken out to ten places 
of decimals. These were obtained in duplicate for each end of the 
quadrant, and the whole were examined by differencing to the sixth 
order. Lastly, the seven-figure logarithms of the functions were 
taken from the eight-figure table of the Service G^ographique de 
l’Arm^e (Paris, 1891), reference being made to Vega where the 
eighth figure was approximately five. 

The log Minus F has been increased by 10 as customary, to 
avoid the inconvenience of printing negative characteristics. 

The whole has been examined by Mr J. Abner Sprigge, of H.M. 
Nautical Almanac Office, so as to make it possible to use the 
tables with confidence in their accuracy to the seventh place. 


(Iota). 

Log plus F. 

A! 

• 

I 

o’0000023 

+ 69 

2 

0092 

115 

3 

0207 

160 

4 

0367 

207 

5 

0574 

253 

6 

0827 

299 

7 

1126 

344 

8 

0*0001470 

+ 390 


a 2 

+ 46 

45 
47 

46 
46 

45 
+ 46 


Log minus F. 
9-9999968 
9871 

9711 

9486 

9196 

8842 

8424 

9*9997942 


Ai 


- 97 

160 
225 
290 

354 

418 

482 

-547 


a 2 


-63 

65 

65 

64 

64 

64 


-65 
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